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Abstract 

Hybrid mesons are exotic mesons in which the color field is not in the ground state. Their 
, understanding deserves interest from a theoretical point of view, because it is intimately related to 

o 

CN ' nonperturbative aspects of QCD. Moreover, it seems that some recently detected particles, such 

p |. 

• as the -7ri(1600) and the y(4260), are serious hybrid candidates. In this work, we investigate the 



(N 



description of such exotic hadrons by applying the auxiliary fields technique (also known as the 
einbein field method) to the widely used spinless Salpeter Hamiltonian with appropriate linear 
^ ■ confinement. Instead of the usual numerical resolution, this technique allows to find simplified 

m ■ 

00 . analytical mass spectra and wave functions of the Hamiltonian, which still lead to reliable quali- 

o , 

I tative predictions. We analyse and compare two different descriptions of hybrid mesons, namely 

O 

^ ' a two-body qq system with an excited flux tube, or a three-body qqg system. We also compute 

r~| ' the masses of the 1 ^ hybrids. Our results are shown to be in satisfactory agreement with lattice 



Oh! 



X 



QCD and other effective models. 
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I. INTRODUCTION 



The study of hybrid mesons is an active domain in theoretical as well as in experimental 
particle physics. From a theoretical point of view, these particles are interpreted as mesons 
in which the color field is in an excited state. Clearly, this problem is related to fundamental 
aspects of QCD, such as its nonperturbative nature. Numerous lattice QCD calculations 
have been devoted to the study of hybrid mesons, in particular to the properties of the l""*" 
state, which is the lightest hybrid with exotic quantum numbers (see Ref. Q for a review. 



and Ref. |2] for more recent_^references) . On ;he experimental side, we can mention the 
recently observed 7ri(1600) 3] and y(4260) which could be interpreted as a 1 ^ nn 
hybrid and a 1 cc hybrid respectively P]. 

Apart from lattice QCD, hybrid mesons have been studied with effective models for a long 
time. For example, we can quote the flux tube model l6|, models with constituent gluons 

Q n 

[21], or the MIT bag model j8||. In potential approaches, to which our paper is devoted, there 
are two main models. In the flrst one, the quark and the antiquark are linked by a string, 
or flux tube, which simulates the exchange of gluons responsible for the confinement. If the 
string is in the ground state, it reduces to the usual linear confinement potential for heavy 
quarks, and to a more general fiux tube model for light quarks, where the dynamics of the 
string cannot be neglected 0, Q|. In this stringy picture, it is possible for the fiux tube 
to fiuctuate, and thus to be in an excited state. These string excitations are analog to the 
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gluon field excitations in full QCD. They have been studied for example in Refs. 
The second approach is to suppose that the hybrid meson is a three-body system, formed 
of a quark, an antiquark, and a constituent gluon, which represents the gluonic excitation. 
Two fundamental strings then link the gluon to the quark and to the antiquark. This picture 
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has been studied in Ref. 13], but also in more recent works 

Nowadays, the spinless Salpeter Hamiltonian (SSH) with a linear confinement is a widely 
used and successful framework to compute hadron spectra in potential models (see previous 
references). Since its kinetic operator is semi-relativistic, most of the calculations have to 
be numerical. However, the auxiliary field (AF) technique allows to greatly simplify the 
calculations [18] and, as we will see, to find analytical solutions to this problem. Even if 
they are approximated, they lead to conclusions which are qualitatively in agreement with 
well-known experimental facts. In particular, Regge trajectories are easily obtained for light 
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mesons by using AF Our purpose is to apply here this formahsm in order to get 

informations about hybrid mesons. This formahsm, also known as the einbein field method, 
has been applied to quark-antiquark two-body systems in Ref. [2^ 21 1. This method can 
be generalized for the case of spinning particles j22 1 . 

Our paper is organized as follows. In Sec. m we solve the SSH in the case of a two-body 
problem. Although this simple case is relatively well-known, it will allow us to introduce 
the AF formalism, and to observe that it leads to correct predictions. Then we present in 
Sec. mil the description of a hybrid meson in terms of a qq system in which the flux tube 
is excited. An other possible approach is to see the hybrid mesons as a three-body system 
made of a quark-antiquark pair and a constituent gluon. This case is studied in Sec. IIVI 
As a result, we are able to compute the effective qq potential in both approaches. These 
potentials can be compared to the predictions of other effective models and of lattice QCD 
calculations. It is done in Sec. |3 Finally, we compute the spectrum of the 1 ^ hybrids in 
Sec. IVH and we sum up our results in Sec. IVIII 

II. THE TWO-BODY PROBLEM 
A. Mass formula and wave function 

The SSH for a system made of two hadrons interacting through a linear confinement is 
given by 



H = ypf + ml + yp^ + ml + ar. (1) 

Let us now introduce three AF (or einbein fields): Two for the quarks, denoted /x,, and one 
for the potential, u. Hamiltonian then becomes 

= ^^ + ^ + ^i±^ + ^ + ^ + ^. (2) 

' ^ 2/ii 2 2^2 2 2z/ 2 ^ ^ 

The AF were introduced to get rid of the square roots in H. Although being formally 
simpler, H{fii, v) is equivalent to H up to the elimination of the auxiliary fields thanks 
to the constraints 



6f,M{fii,u) = ^ /iio = ypf + (3a) 
5yH{jii, u) = Q ^ uq = ar. (3b) 
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It is worth mention ing^ that (/ijo) can be seen as a dynamical mass of the quark whose 
current mass is [23]. Moreover, the Hamiltonian Q can be compared to the one of 
the rotating string model (RSM) [2J]. This is an effective meson model derived from the 
QCD Lagrangian, in which the quark and the antiquark are linked by a Nambu-Goto string 
simulating the confining interaction. The RSM Hamiltonian reads, in the center of mass 
frame, 

^ [ fJ'i Jo Jo aar^J 

with 

as = - 0' + f^2e + [ dPiP- V. (5) 

Jo 

The parameter /3 labels the points of the string. C, defines the position of the center of 
mass on the string, and yUj, v are the AF. In this framework, v is seen as the dynamical 
energy of the string whose "static" energy is ar. The complexity of (j3)) is due to the fact 
that the string contributes to the total angular momentum L. If we neglect the dynamical 
effects of the string, which are in fact sufficiently small to be be treated in perturbation (see 
Ref Q]), Hamiltonian (j^ becomes 

if(/..,.) = ^^ + ^ + ^^ + ^ + ^ + ^, (6) 
' 2/ii 2 2/i2 2 2z/ 2' ^ ' 

which is precisely our SSH Q. 

We can observe from the relations (jH)) that the AF are, strictly speaking, operators. 

However, the calculations are considerably simplified if one considers them as real numbers. 

The elimination of the AF is then finally achieved by minimising the masses with respect 

to them [isj. This procedure leads to a spectrum which is an upper bound of the "true 



spectrum" (computed without AF), the differences being about 10% j25|. In our case, the 
SSH turns out to be a simple nonrelativistic harmonic oscillator (jH]). Its mass spectrum and 
wave functions are thus readily computed. They read 

Af(ft. .) = .(2„ + I + 3/2) + |i + |i + ^^i±|i±i:. (7) 

^ = 0„,,(r)rr(^,^), (8) 

with 



10 = \/ a? / liiy, f3 = \/ fia'^/ u, fi = ^^^^ . (9) 
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Y^™ are the spherical harmonics, and 

= /^^e^+^Z^) v/2n!/r(n + £ + 3/2) e'^'-'/^ L^^ (/Jr^) (10) 

is a properly normalised radial eigenfunction of the three dimensional harmonic oscilla- 
tor j2a]. L° are the Laguerre polynomials. 

V is eliminated by demanding 5yM{jj,i, u) = 0, which leads to 

zyo = (aV^)i/3(2n + ^ + 3/2)2/3^ (11) 

The remaining AF, /ij, cannot be analytically eliminated in general from the condition 
5^^M(/ij, z/q) = 0. We will only consider three relevant special cases. Firstly, if the two 
bodies have a large mass, we can set /ij = because the dynamical effects will be very 
small, and we obtain 

M,, = -^ (2n + £ + 3/2)2/3 + mi +m2, (13) 

with m the reduced mass. The mass formula (fTSj) is valid for example for a meson formed 
of two heavy quarks. It is equal at large £ to the corresponding classical solution of the 
relativistic flux tube model [lO| . This phenomenological model is in fact classically equivalent 
to the RSM if the auxiliary fields are properly eliminated 2^. Moreover, at large angular 
momentum, Mhh oc i"^^^ in qualitative agreement with the experimental data Q] . Secondly, 
if rrii = 0, as it is the case for light mesons and glueballs formed of two gluons, we can 
compute that 

/^i,o = /i2,o = /io = (^) {2n + £ + 3/2)1/2, (^4) 

Ma = 4 (0 (2n + £ + 3/2)^^ = 4/io, (15) 
as it is expected from the relativistic virial theorem j^^]. Squaring p5|l . we get 

= 8a(2n + £ + 3/2). (16) 



When £ is large, it appears that the square mass increases linearly with i. Thus, our solution 
reproduces the Regge trajectories, which are the best known experimental fact concerning 
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the light meson spectroscopy. The Regge slope is given by 8a, which is in agreement with 



a recent calculation of the glueball spectrum with the RSM 



However, it is larger than 
This is in fact related to the 



the prediction of the relativistic flux tube, that is 27ra 
AF technique itself, and more precisely to the number of AF which have to be introduced, 
as explained in Appendix El With a ~ 0.2 GeV^, a mass formula such as (fT!)|l is able to 
correctly reproduce the experimental Regge slope of the mesons P|. Finally if, say, mi = 
and m2 is large, we can find 

(M/,z-m2)' = 4a(2n + £ + 3/2). (17) 



The Regge slope for a meson formed of a li 



two light quarks, as it was shown is Ref. 



iignt 



ht and a heavy quark is thus the half of the one for 



in agreement with experimental observations. 



III. HYBRID MESON AND THE EXCITED FLUX TUBE 

If the color field is in the ground state, it is generally accepted that the potential between 
the quark and the antiquark in a meson is mainly compatible with a funnel potential, 

V,,{r) = ar-^, (18) 

where as is the strong coupling constant. The ar part is pure flux tube, thus pure con- 
finement, while the Coulomb term comes from the one gluon exchange process (OGE). The 
spectrum obtained with (fT^ is in good agreement with experimental data for the light 



and heavy mesons [23], but also with lattice QCD calculations |30j. Typical values for the 
parameters fitting these lattice QCD data are a = 0.2 GeV^, and as = 0.2 — 0.3. 

In a hybrid meson, we have to wonder about how this potential will be modified. A well- 
known approach, based on the computation of the flux tube fluctuations at the quantum 
level, leads to the so-called Liischer term. In this approach, the potential between two fixed 
quarks is given by [3] 

V,^ir) = «^ + 7 - ^) . (19) 

where is the excitation number of the string. For A^ = 0, we recover the Funnel poten- 
tial ()18j) with formally as = n/lG ~ 0.2. This corresponds more or less to the usual value. 
For A^ > 0, the short-range term becomes repulsive, and this potential should become ap- 
plicable to heavy hybrid mesons. 



Another formalism has also been developed in Ref. [ll| to treat the excitations of the 
flux tube. It is based on less conventional approaches to string theory, and leads to 

Vgg{r) = Va^r'' + 27iaN + ^, (20) 

or 

where the short-range term is not due to the string, but again to OGE, with the qq pair in 
an octet. Let us note that for large r, \J a^r'^ + 2TTaN ^ ar + irN/r, and we approximately 
recover the Liischer term. As we can consider that the short-range term can be added 
in perturbation in a first approach, the unperturbed spectrum of the potential ()19j) will be 
defined by the mass formula we derived in Sec. |n]from Eq. (fT^ . However, the confining part 
of potential ()20|) being different, it will affect the mass spectrum, even at the unperturbed 
level. If 7^ 0, the calculations cannot be analytically performed, except for two heavy 
quarks. In this case, one can readily obtain a sort of Regge trajectory with respect to A^, 

(Mhh - Imf ^ 27raA^. (21) 



IV. HYBRID MESONS WITH CONSTITUENT GLUONS 

In this picture, it is assumed that the excitations of the gluon field can be described by the 
potential created by a constituent gluon. The quark- ant iquark pair is thus in a color octet 
in order for the hybrid to be a colorless object. Assuming the Casimir scaling hypothesis, 
which seems to be confirmed by several models js^, it can be shown that the confinement 
is no more a Y-junction like in a baryon but two fundamental strings linking each quarks to 
the gluon [l^. Neglecting all the short-range interactions, the three-body SSH is thus 

i^o = ^ \lvl + m\ + o\xg - Xq\ + a\xg - Xq\. (22) 

Let us now consider a hybrid meson in which the quark and the antiquark have the same 
mass, and are assumed to be static like in lattice QCD. Then, all the properties of the hybrid 
meson should depend only on the quark-antiquark separation and on the quantum numbers 
of the gluon. If we define Xq — Xq = R, and r = —R/2 + Xg — Xq = R/2 + Xg — Xg, the 
potential of the strings is proportional to \r + R/2\ + |r — R/2\. As this expression is not 
useful in a practical computation, we will use the Born-Oppenheimer approximation 

\r + R/2\ + \r- R/2\ ^ 2^ r'^ + i?V4, (23) 
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valid for cos{r,R) <^r/R + R/Ar. This upper bound is always greater than 1 (the mimi- 
mum being reached for 2r = R). Actually, our assumption says that the gluon lies on the 
symmetry plane of the qq pair. The Hamiltonian of the system for static quarks [p^ ^ = 0) 
reads then 

Ho = + 2a^/r^ + R^/4: + 2m. (24) 

The string tension a is the same than in the meson case. We assumed an equal mass m for 
the quark and the antiquark and a vanishing current mass for the gluon. 

The eigenvalues of the Hamiltonian ()24j) can be found for any value of the quark- antiquark 
separation by introducing two AF, fi being again the constituent mass of the gluon and u 
the energy of the strings. With these two fields, Eq. ()24|) becomes a harmonic oscillator 
Hamiltonian, 



2a(2n + i + 3/2) ii + u a'^R^ 

Eo /i, iy)-2m = ' ' + ^1— + 26 

Jjlv 2 2v 



Its eigenvalues 



depend on i?, and on the quantum numbers n and I of the gluon. Their eigenfunctions are 
given by (jH)), with 

13 = 2a^J^. (27) 
The constraints S^^uEolfi, u) = lead to 

z/o = A;Vo (28) 



a^i^Vo + '^Vo - = 0, (29) 

with 

k = 2a{2n + i + 3/2). (30) 

Equation (PHj) defines z/q in terms of fiQ, and this last AF is found to be a solution of Eq. ((221) ■ 
This equation can be solved for any value of R thanks to the Cardan method (see Appendix 
B). 

Two limit cases are interesting. The first one is the limit of large R. Then, /i^ ^ k'^/aR, 
V pa ai?, and 

Eo-2m^ aR. (31) 
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The effective current mass for the gluon decreases with the interquark distance. For large 
quark separation, the energy is only given by the flux tube and the potential energy is the 
expected linear confinement. Secondly, if i? = 0, we have /i = z/ = y/k ^ 775 MeV, and 

(Eo-2m)2 = 4(2a)(2n + £ + 3/2). (32) 

By comparing this formula with (fTTjl . we see that this situation corresponds to a gluelump: 
A hybrid meson seen as a bound state of a gluon and a pointlike heavy meson. The string 
tension is 2a because it is the superposition of two fundamental strings. It can be computed 
that 

(,.^> = (!!i±i±M. (33) 

Consequently, thanks to definition (f^fj) . we have 

^,^.ij,±i±mri^ (34) 



and formula (jH2|) is thus equivalent to 



^0 - 2m = 2 (2a) • (35) 

Half of the energy is given by the confinement and the other half by the kinetic energy of 
the gluon in agreement with the virial theorem. It also very interesting to compare formulas 
(PT|) and (jn21)- They both predict Regge trajectories depending on the color field excitation, 
but they differ in the interpretation they give to it. In the two-body case, this excitation 
is characterised by the quantum number N, defining the state of the string, while in the 
constituent gluon model, the excitation is represented by the gluon itself, and thus the Regge 
trajectory depends on its quantum numbers. 

In the general case, an analytical solution of Eq. ()29|1 can always be found, as shown in 
formula ()B5jl of Appendix El . It is worth noting that our procedure is not the same as 
the one developed in Ref. (l2|, where the AF were eliminated before computing the mass 
spectrum, and the resulting Hamiltonian was solved variationally with a numerical method. 

The confinement interaction gives the correct behavior a large R. But, in order to be 
consistent in the region R <^ ^/a we must add a OGE interaction between each pair of 
particles, i. e. 

AH = 2 ^ ^ + — ^, (36 



with K = —3/2 the color factor of (anti)quark-gluon pair, and k' = 1/6 the color factor 
for the quark-antiquark pair ^l- We used in AH the approximation (j2Sl); as for the unper- 
turbed Hamiltonian (j22I)- The angular momentum is here a good quantum number even 
if we add the short-range interaction. This was not the case in |l2}]. As AH is computed 
perturbatively, the total energy reads, thanks to an usual approximation, 

E-2m^ Eoifio, ^o) - 2m - + ^- (37) 

(r^) is given by Eq. ()34|1 . Let us notice that the Coulomb interaction for the qq system in a 
octet is repulsive. 



V. THE EFFECTIVE QUARK-ANTIQUARK POTENTIAL 

One of the observables in lattice QCD is the potential energy between the static quark- 
antiquark pair. It appears that there are several levels of potential energy, corresponding to 
different states of the gluon field 3^. These excited states of the gluonic field are labeled by 
three quantum numbers. The first one is the projection A of the total angular momentum 
Jg = Lg + Sg of thc gluou ou the qq axis. The capital greek letters S, 11, A, . . . are used to 
indicate the states of A = 0, 1, 2, . . . respectively. The combined operations of the charge 
conjugation and the spatial inversion of the quark and of the antiquark is also a symmetry. 
Its eigenvalue is denoted by rjcp- States with rjcp = 1(— 1) are denoted by the subscripts 
g {u). There is a additional label for the S states: S states which are even (odd) under a 
reflexion in a plane containing the qq axis are denoted by a superscript + (— ). All these 
different states have been observed in Ref. 

In the excited flux tube picture, the glue states with = and A^ = 1 are uniquely 
the ground state E+ and n„ respectively. For A^ > 1, the flux tube can be excited in 



A = 0, 1, . . . , A^ states 



III, 



ll|. The CP value is given by rjcp = (—1)^- Potentials (fT^ and 
are compared to the lattice data in Fig. ^ for the lowest states. As remarked in Ref. 
the particular string potential (j^n|) fit the excited levels with a good accuracy, while the 
Liischer potential (fTTHl reproduces very well the ground state but diverges too fast at small 
i? for A^ > 0. 

In our constituent gluon model, defined by the mass formula (jSZj), we will assume for 
the gluon that n = 0. i is then the only relevant quantum number. Since the qq system 



10 



Luscher 



-4 - 



0.5 



1.5 



2.5 



3.5 



4.5 



5.5 



FIG. 1: Comparison between lattice QCD calculations (symbols) from Ref. [3j|, and two-body 
models like the excited flux tube (solid lines) from Eq. (|2n|) and the Liischer term (dot-dashed 
lines) from (|19|) . All the potentials are plotted in terms of the lattice scale Rq = 2.5 GeV^^ and 
are shifted by an overall amount Vqq{2Ro). The parameters a = 0.2 GeV^ and a = 0.3 are fitted 



and the gluon have both a negative intrinsic parity, the parity of the states is the space 
parity (—1)^. The charge conjugation on the gluon give a —1 factor, and a (— 1)'^«« factor 
for the quark-antiquark pair. The value of rjcp is thus ( — l)^+'599+i and can give either g or 
u states. But, in our case, we made a strong symmetry assumption, and considered that 
the gluon was always located in the symmetry plane of the qq pair. Thus, the only value 
we can reproduce is A = with a positive subscript, since it is the only value in agreement 
with the symmetry of our wave function. With these considerations, we can conclude that 
the states described by our model are S+ states. The spin effect being weaker for maximum 
spin states we should better reproduce the states with Sgg = 1. More precisely, we get 
the S+ for the S-wave {£ = 0), and the for the P-wave {£ = 1). The comparison between 
the lattice S states and the model for i = 0,1 with all short-range interactions is shown in 
Fig. El Although our formula is simple and approximated, one can see that it fits with a 
good agreement the lattice data. Let us note that, following our discussion on the quantum 
numbers C and P, the P-wave corresponds to a l""*" hybrid meson. 

Other recent studies of hybrid mesons with a constituent gluon can also be mentioned 
H H a. In an t.e. wo*, based o„ QCD ,n .Ke Coulo.b .au.e t.e 




on the lattice ground state S^. 
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FIG. 2: Same as Fig. ^ but the solid lines now represent our qqg mass formula (|37|) . The upper 
curve is the P-wave {i = 1), and the lower curve is the S-wave {£ = 0). The ground state is shown 
only for comparison, and is a funnel potential with the parameters of Fig. ^ 




FIG. 3: Comparison between the constituent gluon model of Ref. (symbols), and our qqg model 
without short range interaction (solid lines), with qq OGE term (dot-dashed lines) and with the 
total OGE term (dashed lines). Parameters are those of Fig. ^ 



particle representation, the Hamiltonian of the system was solved variationally by a numer- 
ical method. In a first work 3^ , the level ordering of the states at short distance disagreed 
with lattice results. For example, the Un state was found below the Uu at short distance. Re- 

n 

cently, the authors added new interactions, listed in the diagrams below [1^. a and b are the 
self-energies of the quarks and the gluon respectively, and b and c are the quark-antiquark, 
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a)- 



bj- 



d) 



and (anti)quark-gluon interactions, which are responsible for both the confinement and the 
short-range term. These four diagrams were the only ones considered in Ref. js^. The new 
diagrams are interpreted as particular three-body interactions, and are represented by e and 



/ 



IJ]. We present in Fig. IHlthe comparison between our analytic model and the numerical 



data of Ref. IJ]. The solid lines represent the energy ()2(i|l of the Hamiltonian with only the 
confinement for £ = 0, 1. The dashed-dotted and dashed lines are the energies with the qq 
Coulomb repulsive interaction and with all the short-range interactions respectively. Our 
results are somewhat different at small i?, but agree in the confining part. It can also be 
checked that the full energy (dashed lines) better fits the plot of the first study j^]. This 
can be understood by noting that our model only includes the confinement and the coulomb 
term coming from OGE, thus the diagrams h and c. We neglected the self-energy a and fe, 
but it should not be a dominant contribution in potential models. Indeed, in the case of 



heavy quarks, it is very small 
gluonic self-energy is vanishing 



IS 



Moreover it is even argued in some approaches that the 

nn 

|35l. I36| . However, it seems that the three-body interactions 



e and /, that we did not took into account, bring a relevant contribution to the total energy. 
This could be clarified in future works. 



VI. HYBRID SPECTRUM 

Masses of heavy hybrid mesons can be derived from Eq. ()37|1 by taking the miminum of 
the energy respect to R. We find a that E — 1m = 1.700 GeV is the minimum value for 
the P-wave. As we argued in the previous section, this P-wave represents a 1 hybrid. 
Many lattice studies have been devoted to such exotic mesons, and we can compare our 
predictions with these results. This is done in Tab. IH For the minimum of the S-wave, we 
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find the value E — 2m = 1.16 GeV. This gives a ccg mass around 3.66 GeV for the 1+"*" 
state in disagreement with the lattice result 7.315 GeV of ref. js^]- Nevertheless, our value 
is close to the experimal state X(3872) |38| which could be interpreted as a hybrid meson 
but also as a meson-antimeson molecule 13911. 



Q 


nig 


Mm 


Miat 


n 


0.005 ± 0.003 


1.710 ±0.006 


1.740 ± 0.240 [40] 


s 


0.105 ±0.025 


1.910 ±0.050 


2.100 ±0.120 [41j 


c 


1.250 ±0.100 


4.200 ± 0.200 


4.405 ± 0.038 [2] 


b 


4.500 ± 0.400 


10.700 ± 0.800 


10.977 ±0.123 [42] 



TABLE I: 1 ^ hybrids masses in our constituent gluon model, Mgqg, compared to lattice QCD 
computations, Miat- Masses are in GeV, and the quark masses niq are taken from the Particle 
Data Group 4^1- The errors on our results are simply computed from the Particle Data Group 
errors on the quark masses. 



Obviously, our mass formula is only valid for the heavy quarks, namely c and b. But 
we see in Tab. H] that these masses are even in good agreement with recent lattice results 
concerning light hybrids nng and ssg. This could be an indication that the dominant degree 
of freedom in a hybrid meson is the one coming from the constituent gluon, the quarks 
mainly bringing a rest mass term to the total mass. However, if we compare our results 
to a previous work based on QCD in Coulomb gauge j^l, we find that the agreement is 
good for the charmed hybrid, where the dynamical effects are small, but rather bad in the 
light sector. This shows that the dynamical effects of the quarks, following the way they 
are taken into account, considerably affect the hybrid spectrum. Such effects should thus be 
studied with more details in future works. Let us finally remark that, in our model, we found 
E — 2m = 1.700 GeV to be a constant, which gives the energy coming from the constituent 



gluon in the P-wave. This is close to the exotic meson 7ri(1600) 



In the excited flux 



tube model, the equivalent energy of the string fluctuation is V^naN, with = 3, that is 
1.942 GeV, a higher value than with our constituent gluon approach. In the Coulomb gauge 
model, E — 2m = 1.825 GeV in the charmed sector, in agreement with our result. 
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VII. SUMMARY OF THE RESULTS 



In this work, we studied the hybrid mesons in the framework of potential models. In par- 
ticular, we applied the auxiliary fields technique to obtain analytical mass formula and wave 
functions of the spinless Salpeter Hamiltonian with linear confinement. Firstly, we showed 
that in the well-known case of mesons, the auxiliary fields allow to get easily mass formula 
whose features are qualitatively in agreement with experiment: The Regge trajectories for 
mesons formed of two light quarks or one heavy and one light quark are correctly predicted. 

The simplest way to study hybrid mesons is to work in the excited flux tube framework. 
It is based on the idea that the flux tube (a Nambu-Goto string) linking the quark and the 
antiquark is not in its ground state, but in an excited one. Since the work of Isgur and 

were widely discussed in the litterature. 
In particular, it was shown in Ref. 12] that the first order correction of the excited flux tube 
was a universal term given by 7r(A^ — l/12)/r. In Ref. Q], it was suggested from string 
theory that the confinement potential should be modified to \/ a?r'^ + 2TTaN, this formula 
giving good results when compared to lattice calculations 32]. As in this approach we 
always deal with a two-body problem, the mass spectra can again be easily computed. In 
particular, we showed that in a hybrid meson formed of two heavy quarks of mass m, one 
should have {Mhh — 2m)^ oc A^. 

We also considered an other picture, which assimilates the hybrid meson to a three- 
body quark-antiquark-constituent gluon bound state In the case of two static quarks 
of the same mass, we computed the mass spectrum of the corresponding hybrid meson. 
As a result, we have been able to find an analytic expression for the interquark potential 
in terms of the quantum numbers of the gluon. Although our model was very simple, it 
correctly leads to the gluelump spectrum if the quark-antiquark separation is zero, and 
it reproduces rather well the lattice data which can be described within our assumptions. 
We was actually only able to reproduce the S+ curves, but the other states can also be 
described following the quantum numbers of the gluon (l2| . In the gluelump case moreover, 
one can observe that {Mhh — ^m) oc i. This illustrates the similarity between the excited flux 
tube and the constituent gluon formalisms, but the degrees of freedom are different: The 
excitation number is replaced by the gluon orbital momentum i. Other works studied 



the picture of a constituent gluon with two static quarks [IJ, 12^. An interesting point 
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is that our curves correctly fit previous works 3J], but slightly disagree with more recent 
references, where three-body interactions, which we neglected, are taken into account jl^ . 
The relevance of such interactions in potential models could be studied in the future. We 
have also shown that our model could well reproduce lattice QCD data, and in particular 
the P-wave, corresponding to the 1""*" hybrids. The hybrid masses we obtained are roughly 
in agreement with lattice, and also compatible with the 7ri(1600). However, the mass we 
get for the nn hybrid disagrees with a previous study based on QCD in the Coulomb gauge, 
because of the dynamical effects of the quarks, which are rather strong in this last approach. 

Our study only aims to understand in an intuitive way qualitative features of hybrid 
mesons. In future works, the main challenge will be to compute precisely the various inter- 
actions coming from the dynamics and spin of the particles, which we neglected here. To 
do this, we think that a model with constituent gluons is more interesting, because a con- 
stituent gluon has its well defined color, spin and interactions with quarks. On the contrary, 
the only clear characteristic of an excited flux tube is the quantum number N . We leave 
such detailed studies of the spin effects in hybrids for future work. 

The authors would like to thank the FNRS Belgium and IISN for financial support. We 
are grateful to Dr. Fabian Brau and Dr. Claude Semay for advices and useful discussions, 
and to Dr. Adam. Szczepaniak for providing us the data of Refs. iJ. Is^. 

APPENDIX A: INFLUENCE OF THE AUXILIARY FIELDS ON THE MASS 
SPECTRUM 

Let us suppose that the Hamiltonian of our problem is of the form 

N 

iJo = $^^., (Al) 

i=l 

where Ai are some operators. The eigenvalues Eq and eigenstates \%Ijq) are assumed to be 
known. We can introduce k AF, denoted as 0j, to obtain a new Hamiltonian if^, 

1=1 ^ ^ j=k+l 

As an illustration, one can compare Hamiltonians (jAlj) and ()A2j) to the SSH (0) and (jSj). 
Hk is equivalent to Hq if the AF are directly eliminated as operators. However, our method 
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is to consider that the 0j are real numbers, in order to simphfy the calculations. Let us 
suppose that we know and {ipk) the solutions of the eigenequation defined by H^. Then, 
clearly. 

This implies that 

Ek = {^k\ Hk \i^k) > (^fcl Ho 1^,) > Eo. (A4) 
Moreover, the same argument immediately allows to show that 

Ek+i > Ek- — (^fc+il i^k+i - 4>k+if l^fc+i) 
20fc+i 

> Ek, (A5) 

the only condition being that the (pj are positive. This is always the case in the situations we 
treated here, since the AF are interpreted as effective quark mass and string energy. Thus, 
we have proved that the masses obtained with the AF technique are an upper bound of the 
true masses, as already shown in Ref. 13], but also that this upper bound is less and less 
strong when the number of AF increases. As an example, we can note that the Regge slope 
coming from the full two-body SSH without AF is around 2'ira With only one AF, it 
has been shown in Ref. jisl] that the slope was roughly equal to 7a. Here, with two AF, we 
observed a slope given by 8a. So, the slope increases with the number of AF. Instead of the 
mass formula we should actually write the following inequality: 

Mfi<8a{2n + i + 3/2). (A6) 

APPENDIX B: CARDAN METHOD FOR THIRD DEGREE EQUATIONS 

In this section, we give the explicit solution of Eq. With fi^ = X, this condition can 
be rewritten as 

Since /iq is interpreted as an effective quark mass, X has to be a positive real number. The 
corresponding solution of Eq. (jBljl is 



^0 = V ^ + VQ^ + V^+yV- + V^, (B2) 
17 



with 



Then, we have simply 



« = -SWP (B3) 



2a2i?2 27 V a2/?2 



/^o = +v^. (B5) 



[1] C. Mc Neile, Nucl. Phys. A 711, (2002) 303 |hep-lat/ 020700r, and references therein. 
[2] Y. Liu and X.-Q. Luo, Phys. Rev. D 73, (2006) 054510 |hep-lat/0511015| , and references 
therein. 

[3] G. S. Adams et al, Phys. Rev. Lett. 81, (1998) 5760; S. U. Chung et al, Phys. Rev. D 65, 
(2002) 072001. 

[4] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 95, (2005) 142001 |hep-e x/0506"08r. 
[5] F.E. Close and P. R. Page, Phys. Lett. B 628, (2005) 215 |hep-ph/0507T99l . 
[6] N. Isgur and J. Baton, Phys. Lett. B 124,(1983) 247 ; F. E. Close and P. R. Page, Nucl. Phys. 
B 433, (1995) 233. 

[7] D. Horn and J. Mandula, Phys. Rev. D 17, (1978) 898; A. Le Yaouanc, L. Oliver, O. Pene, 
J.-C. Raynal and S. Ono, Z. Phys. C 28, (1985) 309; Yu. S. Kalashnikova, Z. Phys. C 62, 
(1994) 323. 

[8] T. Barnes and F. E. Close, Phys. Lett. B 116, (1982) 365. 
[9] D. LaCourse and M. G. Olsson, Phys. Rev. D 39, (1989) 2751. 
[10] M. G. Olsson and S. Veseh, Phys. Rev. D 51, (1995) 3578. 

[11] T. J. Allen, M. G. Olsson, and S. Veseh, Phys. Lett. B 434, (1998) 110 |hep-ph/9804"452l. 
[12] M. Liischer, K. Symanzik, and P. Weisz, Nucl. Phys. B 173, (1980) 365; M. Liischer, and P. 

Weisz, JHEP 0207, (2002) 049. 
[13] A. P. Szczepaniak and P. Krupinski, Phys. Rev. D 73, (2006) 034022 |hep-ph/0511083l . 
[14] A. P. Szczepaniak and P. Krupinski, Phys. Rev. D 73, (2006) 116002 |hep-ph/0604098l . 
[15] V. Mathieu, C. Semay and F. Brau, Eur. Phys. J. A 27, (2006) 225 jhep-ph/05112101 . 
[16] E. Abreu and P. Bicudo, pIep^h/050828T| . 



18 



[17] Y. S. Kalashnikova and D. S. Kuzmenko, Phys. Atom. Nucl. 66, (2003) 955 (Yad. Fiz. 66, 
(2003) 988) ||hep-ph/0203" T28l; 67, (2004) 538 (Yad. Fiz. 67, (2004) 556) hep-ph/0302070, ; 
AIP Conf. Proc. 619, (2002) 241 |hep-ph/0110T85] ; Y. S. Kalashnikova and Y. B. Yufryakov, 
Phys. Atom. Nucl. 60, (1997) 307 (Yad. Fiz. 60 N2, (1997) 374) Ihep-ph/gSlOSSTj . 



[18] C. Semay, B. Silvestre-Brac and I. M. Narodetskii, Phys. Rev. D 69, (2004) 014003 
|hep-ph/0309256j. 

[19] F. Buisseret and C. Semay, Phys. Rev. D 71, (2005) 034019 |hep-ph/0412'36n. 



[20] Y. S. Kalashnikova and A. V. Nefediev, Phys. Atom. Nucl. 60, (1997) 1389 (Yad. Fiz. 60, 

(1997) 1529) |hep-ph/961 13611 . 
[21] P. A. M. Dirac, Lectures on Quantum Mechanics, Belter Graduate School of Sciences, Yeshiva 

University, New York (1964) 
[22] L. Brink, P. Di Vecchia and P. S. Howe, Nucl. Phys. B 118, (1977) 76. 



[23] F. Buisseret and C. Semay, Phys. Rev. D 70, (2004) 077501 |hep-ph/0406216] . 

[24] V. L. Morgunov, A. V. Nefediev, and Yu. A. Simonov, Phys. Lett. B 459, (1999) 653. 

[25] W. Lucha and F. F. Schoberl, Phys. Rev. A 54, (1996) 3790 p^ep^h/ 9603429 ; F. Buisseret 

and C. Semay, Phys. Rev. E 71, (2005) 026705 |heF-ph/0409033j . 
[26] F. Brau, J. Phys. A 32, (1999) 7691 |quant-ph/9905033| . 



[27] W. Lucha and F. F. Schoberl, Phys. Rev. Lett. 64, (1990) 2733. 

[28] Yu. A. Simonov, |hep-ph/0603148] . 

[29] S. Godfrey and N. Isgur, Phys. Rev. D 32, (1985) 189. 

[30] Y. Koma, E.M. Ilgenfritz, T. Suzuki, and H. Toki, Phys. Rev. D 64, (2001) 014015 
Ihep-ph/OOlT Tes'; T. T. Takahashi, H. Matsufuru, Y. Nemoto, and H. Suganuma, Phys. 
Rev. Lett. 86, (2001) 18 |hep-lat/0204"oTT| . 

[31] S. Deldar, Phys. Rev. D 62, (2000) 034509 |hep-lat/99H"008] ; G. S. Bah, Phys. Rev. D 62, 
(2000) 114503 ^glat/0006022 ; C. Semay, Eur. Phys. J. A 22, (2004) 353 pIep^ph/0409105]. 

[32] K. J. Juge, J. Kuti and C. J. Morningstar, Nucl. Phys. Proc. Suppl. 63, (1998) 326 



|hep-lat/9709131. ; K. J. Juge, J. Kuti and C. Morningstar, Nucl. Phys. Proc. Suppl. 119, 
(2003) 682 |hep-lat/0209"l09] . 
[33] K. J. Juge, J. Kuti and C. Morningstar, Phys. Rev. Lett. 90, (2003) 161601 |hep-l at/0207"00?. 



[34] E. S. Swanson and A. P. Szczepaniak, Phys. Rev. D 59, (1999) 014035 p^ph/9804219 
[35] Yu. A. Simonov, Phys. Lett. B 515, (2001) 137 [hep- ph/0105141]. 



19 



[36] A.B. Kaidalov, Yu.A. Simonov, Phys. At. Nucl. 63, (2000) 1428 (Yad. Fiz. 63, (2000) 1428) 
|hep-ph/991 1291 1 ■ 

X. Q. Luo and Y. Liu, Phys. Rev. D 74 (2006) 034502 |hep-lat/0512044l . 
S. K. Choi et al. [Belle Collaboration], Phys. Rev. Lett. 91 (2003) 262001 |hep-e x/0309"032|. 
A. A. Petrov, |hep-ph/0 602066 ; L. Maiani, F. Piccinini, A. D. Polosa and V. Riquer, AIP 
Conf. Proc. 814 (2006) 508 hep-ph/05 12082]. 
J. N. Hedditch, Phys. Rev. D 72, (2005) 114507 |hep-lat/0509"l06| . 
C. Bernard et al, Phys. Rev. D 68, (2003) 074505 |hep-lat/0301024] . 
T. Manke et al, Phys. Rev. D 64, (2001) 097505 |hep-lat /01030T5|. 
Particle Data Group, S. Eidelman et al, Phys. Lett. B 592, (2004) 1. 
F. J. Llanes-Estrada and S. R. Cotanch, Phys. Lett. B 504, (2001) 15 hep-ph/0008337, . 
N. Isgur and J. E. Baton, Phys. Rev. D 31, (1985) 2910. 



20 



